Electron Self-Energy and Generalized Drude Formula 
for Infrared Conductivity of Metals 
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Gotze and Wolfle (GW) wrote the conductivity in terms of a memory function A4{uj) as f7{ui) = 
{ine'^ /m){uj+M{uj))~^ , where M{uj) = i/r in the Drude limit. The analytic properties of — M (oj) are 
the same as those of the self-energy of a retarded Green's function. In the approximate treatment of 
GW, —M closely resembles a self-energy, with differences, e.g., the imaginary part is twice too large. 
The correct relation between —M and E is known for the electron-phonon case and is conjectured 
to be similar for other perturbations. When vertex corrections are ignored there is a known relation. 
A derivation using Matsubara temperature Green's functions is given. 



I. PRELIMINARIES 

Holsteini used elementary arguments to show that in 
the infrared properties of metals there can be quan- 
tum effects (outside of the semiclassical Boltzmann ap- 
proach) when the temperature is low enough and the 
probing frequency lo is degenerate with non-electronic 
excitations like phonons. Such effects have been seen 
experimentallySi^. Gotze and Wolfle (GW)^ gave a nice 
simplified way to compute such effects in the optical re- 
sponse of metals using truncated equations of motion to 
compute the "memory function" M{uj) defined as 



cr(w) 



ine^/m 
u + M{w) 



(1) 



In the dc limit, their formulas correctly reproduce lowest- 
order variational solutions of the corresponding Boltz- 
mann transport theories. Unfortunately, a systematic 
perturbation theory for M{uj) is not known. The GW 
results are slightly less accurate than the corresponding 
lowest order results of Green's function theories. 

The function —M{u}) has causal analytic properties 
and, not surprisingly, bears a close resemblance to an 
electron self-energy Yi{k,Lo) for fc-points averaged over 
the Fermi surface. However, the imaginary part of S 
is — 1/2t while the imaginary part of —M must be — l/r. 
This is not the only difference between —M and E. Since 
the analogy between —M and E is sometimes used for 
analysis of infrared spectra^, it is important to under- 
stand just how good it actually is. A full formula seems 
not to have been derived, and is beyond the ambition 
of this paper. A reasonable conjecture is that when 
anisotropy with k around the Fermi surface is not too 
important, then 



a{uj) 



duj'- 



f{uj')-f{0j' + L0) 



muj Lu - Eii.(tj' + Lu + ir/) -|- E* (w' -|- 177) 

(2) 

where /(w') = (exp(/3w') -f 1) ^ is the Fermi-Dirac func- 
tion. Here Ej^ is a modified version of E, averaged over 
the Fermi surface, with a complex weighting factor, simi- 
lar to the familiar transport factor I — cos{9). The actual 
weighting factor (in the solved electron-phonon caseSii) 



is found from a complex, frequency-dependent non-linear 
integral equation. Replacement of the weight factor by I, 
turning E;,. into an ordinary, but A:-averaged, self-energy, 
should work fairly well in most cases. Sher— did a nu- 
merical study which tends to confirm that the difference 
between Eji- and E is small. 

Eq. (O implies a relation between —M{oj) and the 
self-energy which becomes more direct at low frequencies. 
Keeping the lowest order (in to) terms, one gets a deriva- 
tive of the Fermi-Dirac function —df{uj')/duj' which can 
be approximated by the Dirac 5{uj'). 



duj' 



duj' 



LU{1 - dT,ir,l{u)')/duj') + 2iEir,2(w') 

ine^/m 

LU — UjdY.ir,l{L0)/dLU + 2iEir,2(t^) 



(3) 



Therefore the real part of —M at low frequencies is 
Ci;(iReEir(aj)/da;, and the imaginary part is — 2ImEir(ti;). 
If the interesting part of ReE is odd in uj, and approxi- 
mately linear, then M at very low w is a lot like E except 
for the factor of 2 in the imaginary part. 

In the dc limit, the result a = ne'^r/m is retrieved, 
with 1/t — — 2ImEii.(a; — > 0). There are minor differ- 
ences between this and the more exact result found from 
a solution of the Boltzmann transport equation. These 
differences arise from fc-averaging and disappear for a 
spherical Fermi surface. 



II. 



KUBO FORMULA 



The starting point is the Kubc^ formula for the con- 
ductivity. For external electric field E(t) = E cos{Lut) the 
current operator j = —eJ2k '^kxc\cf^ acquires a non-zero 
expectation value j(t) = Re[(T(ti;) exp(— where the 
linear response coefficient 17(0;) is 



r{uj) + 

m 



(4) 



2 



(5) 



The Hamiltonian 7i — TLq + Tii has the non-interactmg 
part Hq — <ikc\cf.. The label k is short for the Bloch 
wavevector and other quantum numbers (kna). The 
state k has energy and group velocity Vk- 

To obtain a Wick-ordered perturbation theory we use 
an imaginary time (0 < a < f3 — l/ksT) version of r{uj), 



dae'^-''{Tj{a)j{0)) 



(6) 



where j(a) — exjp{aTi.)j exp(— crH). Angular brackets de- 
note both a canonical ensemble temperature average and 
if necessary, an average over an ensemble of randomly 
distributed impurities. The Matsubara frequency is 
2nfi/(3 and is an integer. When analytically contin- 
ued from iujfj, to uj + irj just above the real lo axis {rj is a 
positive infinitesimal) r{iuj^) becomes r(a;), the retarded 
correlation function needed for the Kubo formula. 

All Feynman graphs for r{iujf^) are formally summed 
in terms of the exact electron Green's function 



G{k, iuj^) 



1 



- €k- T.{k,iu)^) 



(7) 



and the exact vertex function T{k,k',iLL!f^,iLU,y), where 
— 2-k{v -\- 1/2)/ (3. The exact answer is 



Vk'x'i'ikk',', 



kk' 



X G{k,iuj^ + iuj^)G{k,iuj^) 



(8) 



Neither S nor F can be calculated exactly. A linearized 
Boltzmann equation is obtained when lowest order re- 
sults for E and F are treated consistently. 

An explicit formula relating ct to S occurs when F is 
replaced by its lowest order term, 



T{kk', iuj^, iujjy) ^ Fq = Vkx5{k, k'). 



(9) 



The corresponding answer for cr(w), denoted by (To(ti^), 
after continuing to the real frequency axis, is 



CTo(w) 



dio' - 



Lu — Tj{lu' + lu + irj) + S* (cj' + irj) 

(10) 

This is not an original result, but since derivations are 
obscure, an elementary one is given in the next section. 

Unlike the conjectured version Eq. the approx- 

imation of Eq. ^ does not correctly reproduce the 
Boltzmann dc conductivity because of the omission of 
vertex corrections. The quasiparticle scattering rate 
1/r = — 2ImE differs from the transport scattering rate 
l/^tr = — 2ImEir by the factor "1 — cos 6*." This comes 
from scattering-in terms of the Boltzmann scattering op- 
erator, approximately the same as the contribution from 
the missing "rungs" of the vertex function omitted in 



Fq. The omitted correction takes into account that small 
angle scattering events do not degrade the current ef- 
ficiently and make smaller contributions to 1 /rtr than to 
1/t. The difference, except at low temperatures, is likely 
to be numerically small, since small angle scattering does 
not usually play a dominant role. 



III. DERIVATION OF EQ. (fTHjl 

Starting by inserting Eq. into Eq. © , 



roiii^fi) = '"kxG{k, iuju + ii^^,)G{k, iuj^) (11) 

k 

This approximation, labeled rg, keeps in principle arbi- 
trarily complicated self-energy graphs in G. 
The spectral function is defined as 

G{k,iuj,) = r du^^^ (12) 



A{k, Lo) = ImG(fc, iuji, —^uj + iri) (13) 

TT 

where G{k,uj + irj) is the "retarded" Green's function, 

1 



G(fc, LU + IT]) 



LU - Ck - T.{k,uj) 



(14) 



and S(/c,(x') = A{k,uj) — i/2T{k,Lu) has imaginary part 
1/2t non-negative. The approximate correlation func- 
tion To becomes 

g2 l-oo l-oo 

roii^^ij.) = ~-t2_/^^^ / '^'^ / duj2A{k,u;i)A{k,uj2) 



E 



1 



1 



ia;,y -|- ILU/^i — UJl ILUi, — LU2 

The Matsubara sum can be performed exactly, 

1 1 1 /(^2)-/(^l) 



(15) 



-y 



3 ^ — ' iiUi^ + ILUn — UJl ILUi, — LU9 

The correlation function now is 



0^2 — t^l + ILiJfj^ 

(16) 



/OO nOO 
- OO 7. J — OO 



X 



dLU2A{k,LUl)A{k,LU2) 



/(^2) - fjuJi) 
UJ2- UJl+ iUJn, ' 



(17) 



where a gratuitous factor 1 = J deS{e — Ck) was inserted. 
Because of the S function in the fc-sum 



Y,vLK^ - ^k)A{k,^l)A{k,LU2) 



(18) 



it is allowed to replace Ck in the denominators of the 
spectral functions by e. There is a rapid dependence on 



3 



e in A{k,uj) which must be treated carefully, and a re- 
maining slow k dependence of S(fc, w) in the denominator 
of A{k,oj) which can be treated less carefully. The self- 
energies in the spectral functions — (l/7r)Im(a; — efe — S)^^ 
are replaced by their /c-average over the Fermi surface, 



SM=^S(fc,C.)5(6fe)/^<5(6fe). 



(19) 



This is usually not a bad approximation, because the 
anisotropy of I](fc,ct;), as k varies around the Fermi sur- 
face, is usually small. The errors introduced by doing 
this are of the same nature as the errors introduced by 
neglecting all vertex corrections. The fc-sum is now 



fe 



k 

h' ^ dkl 



f;2 X! 



./ = Km]eff(e).(20) 



Here the b function was replaced by —df /dek- In- 
tegrating by parts gives the inverse effective mass 
[d^tk/ dk'^) /Ti^ summed over all states lower in energy 
than e. 

The range of the remaining e-integration is nominally 
(— cx), oo). However, the factors A(k,oj) are peaked at 
t — LOi and e = L02] loi and uj2 thus have similar values. 
Because of the factor (/(W2) — /(wi)), they must both 
lie near the Fermi energy (one below and one above.) 
Therefore the e integral is dominated by e near the Fermi 
energy. The value of [n/m\ce{e) at the Fermi energy is 
[n/m]eff, the number of electrons divided by the effective 
mass averaged over all states below the Fermi energy. An 
equivalent formula is 



It is necessary to integrate e carefully over the Lorentzian 
peaks of A(e, a;i)A(e, W2)- Cauchy's theorem can be used 
after closing the e-contour by an arc going to infinity 
in either the upper or lower half-plane. The result is 
expressed by another identity, 



rfe (i) Im 



uji - e - Iji J \7r J \uj2 - e - 2^2 



i) Im f I 

Try \LUi - UJ2 - ^1 + ^2 



(23) 



The proof is elementary but tedious. The current corre- 
lation function is now 



' n 
-m 
X Im 



e 

3ff TT 



duji 



°1 



LUi — UJ2 — ^1 + So 



-oox UJ2 - i^l + iuJ^j. 



(24) 



The function ro(w) is now just ro(ia;^) with iw^ re- 
placed by w -|- IT]. The only complex quantity in the for- 
mula for ro(w) is the factor {lu2 — uji + oj + irj)^^, so the 
real part ReCTo(ti') = lnirQ{Lu) /uj is 



e 

off to 



2 roc 



duji I dw2[/(w2) - /(wi)] 
X S{uj2 — ciJi -|- w)Re 



" n ' 












.m. 


off ^ 





duj' 



071 — 0^2 — Si + E; 



X Re 



(25) 



The current correlation function now is 



7'o(«w^) 



de 



duji 



X A(e,c.i)A(e,a;2) ^^"^^"{^"^^ (22) 

OJ2 — OJl + lOJu 



The function ao{Lu) is specified by the requirements of 
being analytic for Imw > 0, vanishing sufhciently rapidly 
as w 00, and agreeing with Eq. (|25|l . It is necessary 
and sufficient to remove the "real part" designator from 
both sides. This is the derivation of Eq. ifTUIl . 
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